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1. Introduction 

The layout design for the Journal of Plasma Physics has been implemented as a 
1^1^25 class file. The JPP class is based on the ARTICLE class as discussed in 
the Wr^K. manual (2nd edition). Commands which differ from the standard I^T5t;X 
interface, or which are provided in addition to the standard interface, are explained 
in this guide. This guide is not a substitute for the I^T^jX manual itself. 

Note that the final printed version of papers will use the Monotype Modern 
typeface rather than the Computer Modern available to authors. Also, the measure 
in JPP class is different from Wr^2s article class. For these reasons line lengths 
and page breaks will change and authors should not insert hard breaks in their text. 

1.1. Introduction to BT^ 

The I^T^X document preparation system is a special version of the TJrjX typesetting 
program. I^TJi?^ adds to Tf^ a collection of commands which allow the author to 
concentrate on the logical structure of the document rather than its visual layout. 

I^TJt;X provides a consistent and comprehensive document preparation interface. 
1^1^^ can automatically number equations, figures, tables, and list entries, as well 
as sections and subsections. Using this numbering system, bibliographic citations, 
page references and cross-references to any other numbered entity (e.g. section, 
equation, figure, list entry) are quite straightforward. 

1.2. General style issues 

Use of defaults will result in a pleasing uniformity of layout and font selection. 

Authors should resist the temptation to make ad hoc changes to these. Also avoid 
use of direct formatting unless really necessary. Papers will be edited as usual, and 
this process may be obstructed by the use of inserted line breaks, etc. 

For general style issues, authors are referred to the 'Preparation of manuscripts' 
in the back cover of the journal. The language used in the journal is British English, 
and spelling should conform to this. 

Use should be made of symbolic references (\ref ) in order to protect against late 
changes of order, etc. 

1.3. Submission of 3TgX articles 

Authors who intend to submit a Mjt;X article should obtain a copy of the JPP class 
file. This is available by anonymous FTP from 

ftp . cup . Ccun .ac.uk 

You will find the class file and instructions contained in a single zip archive jpp-cls . zip 
in the directory 

pub/texarchive/ j ournals/latex/ j pp-cls 

The readme . txt (which is the same directory) tells you how to unzip the jpp-cls . zip 
file. There may also be an 'unpacked' directory containing all of the files separately, 
in case of difficulty. If you cannot obtain the JPP files, use the standard article 
class, with the default 'lOpt' option. 

When submitting the final article, ensure that the following are included and are 
clearly labelled. 

{a) A hardcopy printout of the article. 
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(6) The input file (exactly matching the hardcopy), 
(c) A copy of any user-defined macros. 

(rf) If you have used BibT^jX, the .bib, .bbl and .bst files that were used. 

(e) Any other files necessary to prepare the article for typesetting. 

The files for the final article should be supplied on a PC or Macintosh disk. Please 
write on the disk what format it is, as this saves time. Submit the hardcopy and 
disk by post in the normal way. 

1.^. The JPP document class 

The use of document classes allows a simple change of style (or style option) to 
transform the appearance of your document. The CUP JPP class file preserves the 
standard IMi^X interface such that any document which can be produced using the 
standard I^T^jX ARTICLE class can also be produced with the JPP class. However, 
the measure (or width of text) is slightly different from that for ARTICLE; therefore 
line breaks will change and it is possible that equations may need re-setting. 



2. Using the JPP class file 

First, copy the file jpp.cls and upmath.sty (see §2.2) into the correct subdi- 
rectory on your system. The JPP document style is implemented as a complete 
document class, and not a document class option. In order to use the JPP class, 
replace article by jpp in the Xdocumentclass command at the beginning of your 
document: 

\documentclass{article} 

is replaced by 

\documentclass{jpp} 

Some of the standard document class options should be used. Author-defined mac- 
ros should be inserted before \begin{document}, or in a separate file and should 
be included with the submission, see §5. Authors must not change any of the macro 
definitions or parameters in jpp.cls. 

2.1. Document class options 

In general, the following standard document class options should not be used with 
the JPP class file: 

• lOpt, 11 pt and 12pt - unavailable; 

• twoside is the default (oneside is disabled); 

• onecolumn is the default (twocolumn is disabled); 

• titlepage is not required and is disabled; 

• f leqn and leqno should not be used, and are disabled. 

The standard document class options should not be used with the JPP class: 

• lOpt, llpt, 12pt - unavailable. 

• draft, twoside (no associated class file) - twoside is the default. 

• f leqn, leqno, titlepage, twocolumn - unavailable. 
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2.2. Additional 'packages' supplied with jpp . els 

The following additional package ( . sty) files are supplied in the JPP distribution: 

• upmath - provides 'upright' Greek math characters. Requires the AMS amsbsy 

and amsgen packages. 

If your site does not have the AMS Fonts and AMSM^j?^ packages installed, 
we strongly recommend that your site installs them. With them, you can produce 
output which is much closer to the final result. The latest AMS Fonts/ AMS I^T^X 
distributions can be found on your nearest CTAN (Comprehensive TjjpC Archive 
Network) site. 

2.2.1. The UPMATH package The upmath package defines the macros \mathup and 
Xmathbup, which allow access to the symbols in the AMS Euler fonts. 

The upmath package provides macros for upright lower-case Greek (\ualpha- 
\uxi), and for bold lower-case Greek (\ubalpha-\ubxi). The bold upright symbol 
\eta has to be treated differently, in this case use \uboldeta. The upmath package 
also provides \upartial and \ubpartial. 

To use the upmath package, you need to have the AMS eurm/b fonts installed. 

3. Additional facilities 

In addition to all the standard design elements, the JPP class file includes 

the following features: 

• Extended commands for specifying a short version of the title and author(s) 
for the running headlines. 

• Abstract environment 

• Control of enumerated lists. 

• An extended theorem environment, enabling you to also typeset unnumbered 
theorems, etc. 

• Ability to create new theorem-like environments for claims, conjectures, exam- 
ples, problems, remarks, etc. 

• The subequations and subeqnarray environments for sub-numbering equa- 
tions. 

• A proof environment. 

• A acknowledgements environment. 

Once you have used these additional facilities in your document, it can be processed 
only with jpp. els. 

3.1. Titles, authors' names and running headlines 

In the JPP c;lass, the title of the article and the author's name (or aiithors' names) 
are used both at the beginning of the article for the main title and throughout 
the article as running headlines at the top of every page. The title is used on odd- 
numbered pages (rectos) and the author's name (with initials only for first names) 
appears on even-numbered pages (versos). The \pagestyle and \thispagestyle 
commands should not be used. Similarly, the commands \markright and \markboth 
should not be necessary. 

Although the article title can run to several lines of text, the running headline 
must be a single line. Moreover, the title can also incorporate new line commands 
(e.g. \\) but these are not acceptable in a running headline. To enable you to specify 
an alternative short title and author's name, the standard \title and \author 
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commands have been extended to take an optional argument to be used as the 
running headUne: 

\title[A short title] {The full title which can be as long 

as necessary} 

\author [Initials and last neimes of all authors] {The full names of 
all the authors, using letterspacing and \cind before 
the last name in the list} 

\af f iliation{Author ' s affiliation} 

Unlike most other document classes, j pp. els does not follow the convention 
of using the second line of the \author command to write the author's affilia- 
tion. Rather, this must be entered using the \affiliation command as shown 
above. The \and command does not produce separate author /affiliation pairs, but 
only generates the word 'and' in small caps as required by the JPP class. Use the 
predefined macro \ls to get the letterspacing between the letters of the author's 
name. If there is more than one affiliation, use \\ [\af f ilskip] between lines of 
the \af filiation, inserting footnote numbers manually. 

An example will make all this clearer. To produce 

Interesting title 

ALAN N. OTHERS and J.Q. PUBLIC^ 

^University of Moscow, Moscow, Russia 
■^Ngonga University, Nairobi, Kenya 

(Received 12 June 1992) 

you would type 

\title [Short title] {Interesting title} 
\author[A. N. Other and J. Q. Public] 

{A\ls L\ls A\ls N\ns N.\ns 0\ls T\ls H\ls E\ls R$-l$\ns 
\and\ns J.\,Q.\ns P\ls U\ls B\ls L\ls I\ls C$-2$} 
\af f iliation{$"l$University of Moscow, Moscow, RussiaW [\aff ilskip] 

$"2$Ngonga University, Nairobi, Kenya} 
\date{12 June 1992} 

The M^jX \thanks command can be used inside the \author command to insert 
an authors 'Current address' information as a footnote, appearing at the bottom 
of the title page. 

3.2. Abstract 

The JPP class provides for an abstract, produced by 
\begin{abstract} 

\end{abstract} 

This should appear just before the first \section command. 

3.3. Lists 

The JPP style provides the three standard list environments plus an additional 
unmunbered list: 

• Numbered lists, created using the enumerate environment. 
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• Bulleted lists, created using the itemize environment. 

• Labelled lists, created using the description environment. 

The enumerate environment numbers each list item with an italic letter in paren- 
theses; alternative styles can be achieved by inserting a redefinition of the number 
labelling command after the \begin{enumerate>. For example, a list numbered 
with roman numerals inside parentheses can be produced by the following com- 
mands: 

\begin{ enumerate} 

\renewcommaiid{\theenumi}{\romcin{enumi]-} 
\item first item 

\end{ enumerate} 

This produces the following list: 

(i) first item 

(ii) second item 

(iii) etc. 

In the last example, the item alignment is uneven because the standard list in- 
dention is designed to be sufficient for arable numerals rather than the wider roman 
numerals. In order to enable different labels to be used more easily, the eniimerate 
environment in the JPP style can be given an optional argument which (like a 
standard thebibliography environment) specifies the widest label. For example, 

(i) first item 

(ii) second item 

(iii) etc. 

was produced by the following input: 

\begin{enumerate} [(iii)] 
\renewcommEind{\theenumi}{\roman{enumi}} 
\item first item 

\end{enumerate} 

Once you have used the optional argument on the enumerate environment, do 
not process your document with a standard I^TJjX class file. 

3.4- Theorem environments 

The \newtheorem command works as described in the M^TjX manual, but produces 
spacing and caption typefaces required to the JPP style. The preferred numbering 
scheme is for theorems to be numbered within sections, as 1.1, 1.2, 1.3, etc., but 
other numbering schemes are permissible and may be implemented as described in 
the M^jX manual. In order to allow authors maximum flexibility in numbering and 
naming, no theorem-like environments are defined in jpp. els. Rather, you have to 
define each one yourself. Theorem-like environments include Theorem, Definition, 
Lemma, Corollary, and Proposition. For example, 

\newtheorem{lemma}-CLemma} [section] 
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\renewcoiniiiaiid{\theleinma}-{\Romcin{leinma3-} 

The new macro which follows can also produce numbered and unnumbered en- 
vironments. 

3.5. Definitions 

The \newdef inition command can be used to create environments for claims, 
conjectures, examples, problems, remarks, etc. These are typeset in the way as 
theorems, except the text is typeset in roman instead of italic. 

You can have an environment created by \newdef inition number in the same 
sequence as a 'Theorem' by adding the optional [theorem] argument as below: 

\newdefinition{def inition} [theorem] {Definition} 

\begin{def inition} 

This is a definition. 

\end{def inition} 

\begin{def inition*} 

This is an unnumbered definition. 

\end{def inition*} 

which produces the following: 

Definition 3.1. This is a definition. 

Definition. This is an unnumbered definition. 

3. 6. Proof environment 

The standard constructs do not include a proof environment to follow a 

theorem, lemma etc. (see also §3.4), and so one has been added for the JPP style. 

Note the use in the following examples of an optional argument in square braces 
which may contain any information you may wish to add. For example, 

\begin{theorem} [2] \label{thenv} 

Let the scalar function $T(x,y,t,\biiiath{\omega})$ be a conserved 
density for solutions of {\rm(9)}. Then the two-component function 
\begin-[equation} 

\hbox{\boldit P} = \mathsf -[J}{\cal E} T 
\end{ equation} 

represents the infinitesimal generator of a symmetry group 

for -[\rm (9)}. 

\end{theorem} 

\begin{proof} [Proof of Theorem~\ref {thenv}] 

The assumption about $T$ meams that 

\[ 

\sim \f rac{\partial T}{\partial t} + 
{\cal E}T \bmath{\omega}_t 

= \frac{\partial T}{\partial t} + \{ T, H\}, 
\] 

where $\partial T / \partial t$ refers to explicit 

dependence on $t$. The skew symmetry of $\mathsf{J}$ hence implies 

\begin-[equation} 
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\frac{\partial T}{\partial t> Xsiin \{ H, T\}, 
\end{equation} 

whereupon the operation $\mathsf {J}\cal E$, which commutes with 

$\partial_t$ in its present sense, gives 

\[ 

\frac{\partial \hbox{\boldit P}}{\partial t} 

= \mathsf{J}{\cal E}\{ H, T\}. 

\] 

This equation reproduces the characterisation of symmetries that 
was expressed by (19) , thus showing P to represent a symmetry 
group . 
\end{proof } 

produces the following text: 

Theorem 3.2 (2). Let the scalar function T{x,y,t,u}) be a conserved density 
for solutions of (9). Then the two-component function 

P=]£T (3.1) 

represents the infinitesimal generator of a symmetry group for (9) . 

Proof. [Proof of Theorem 3.2] The assumption about T means that 

dT dT 

where dT/dt refers to explicit dependence on t. The skew symmetry of J hence 
implies 

dT 

-Q^-{H,T}, (3.2) 
whereupon the operation ^£, which commutes with dt in its present sense, gives 

dP 

This equation reproducers the characterisation of symmetries that was expressed by 
(19), thus showing P to represent a symmetry group. □ 

The final □ will not be included if the proof* environment is used. 



4. Mathematics and units 

The JPP class file will insert the correct space above and below displayed maths if 
standard M^rpC commands are used; for example use \ [ ... \] and not $$ ... $$. 
Do not leave blank lines above and below displayed equations unless a new para- 
graph is really intended. 

4-1- Numbering of equations 

The subequations and subeqnarray environments have been incorporated into 

the JPP class file. Using these two environments, you can number your equations 
(1.1a), (1.1&) etc. automatically. For example, you can typeset 

ai = {2nM^/x)^y^ (4.1a) 
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and 



02 = {x/2VL)^ky/M. 



(4.1b) 



by using the subequations environment as follows: 

\begin{subequations} 
\begin{equation} 

a_l \equiv (2\0mega M-2/x) -{\f rac{l}{4}} y-{\f rac{l}{2» 
\end{equation]- 
and 

\begin{equation} 

a_2 \equiv (x/2\0mega) "-CXf rac{l}{2}>k_y/M. 
\end{equation} 
\end{subequations} 

You may also typeset an array such as: 



by using the subeqnarray environment as follows: 

\begin{subeqnarray} 

\dot{Xj k = k \gamma X - \gELinma\delta\eta ,\\ 
\dot{\eta} k = k {\textstyle\f rac{l}-[2}} \delta + 2X\eta . 

\end{ subeqnarray} 

4-2. AMS fonts - especially if you do not have them 

If you need AMS symbols but do not have the AMS fonts you can ensure that 
they will be correctly typeset by taking the following steps. Set up user-defined 
macros that can be redefined by the typesetter to use the correct AMS macros. For 
example, the blackboard bold symbols, sometimes called shell or outline characters, 
are obtained with the AMS macro \mathbb{ . . }. Instead, use a macro definition 
such as: 

y. replace font! 

\newcoinmand\BbbE-[\ensuremath{\mathsf {E}}} % Blackboard bold E 

This substitutes a sans serif character where you want blackboard bold. You can 
typeset the input file and the typesetter is alerted to do the substitution. 

The following example (which uses the Xprovideconmicind macro) will work with- 
out modification by the typesetter, because the Xprovidecommand macro will not 
overwrite any existing \mathbb definition. 

\providecoinincind\mathbb [1] {\ensuremath-C\mathsf {#1}}} 
\newcoinmand\BbbE{\inathbb{E}} Blackboard bold E 

Plain T^TjX provides only \leq and \geq which typeset the Computer Modern 
symbols < and >, respectively. These will be redefined at typesetting to use the 
AMS symbols \geqslant and \leqslant, to give the slanted symbols. 

If you wish to use AMS fonts with WT^^2s you must be using at least version 
2.0. Earlier versions are not supported. 



X = jX — ^St], 



f]=\5 + 2XT^. 



(4.2a) 
(4.2b) 
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4-3. Typeface 

Sometimes, non-italic symbols are required in maths. This section describes how 
these may be obtained using and jpp. els. 

4.3.1. Roman symbols The mathematical operators and constants, such as sin, 
cos, log and cxp (and many others) arc covered by Wr^i^2s macros which ensure 
that they are typeset in roman text, even in math mode: \sin, \cos, \log, \exp. 
Where single letters are concerned (e.g. d, i, e) just use \inathrin in math mode, i.e 
$E=m\mathrm{c}~2$ which typesets as = mc^, giving the correct roman character 
but with math spacing. When the term involves more than one character (e.g. Re 
or Im) text-character spacing is required: 

$\mbox{Re}\;x$ 

which typesets as Re x. 

Where such expressions are used repeatedly, macro definitions can reduce typing 
and editing. The following examples are included in the preamble of the input files 
for this document, jppguide.tex, and of the sample article jppsampl .tex. Authors 

are encouraged to use them and others like them. 

\newcoinincind\Real{\mbox{Re}} °/o do not confuse with TeX's \Re 
\newcoinincind\Iinag{\inbox-[Im}} % do not confuse with TeX's \Iin 

\newcommand\Ai{\mbox{Ai}} 
\newcommand\Bi{\mbox{Bi}} 

4.3.2. Multiletter italic symbols If multiletter symbols are used in maths mode, 
for example Reynolds, Prandtl numbers, etc. the standard maths mode spacing 
between them is too large and text-character spacing is required. As described in 
§4.3.1 (but here for italic letters) use for example 

\newcoininEind\Rey{\mbox{\textit{Re)-}-} 
\newcoinmand\Pran-[\inbox{\textit{Pr}}} 

4.3.3. Upright and sloping Greek symbols Like the italic/roman fonts described 
above, Greek letters for math variables arc printed in the Journal in sloping type, 
whereas constants, operators and units are upright. However, normally pro- 
duces sloping lowercase Greek and upright capitals. The upright lowercase required 
for example for 'pi' and 'mu' (micro) are provided by the upmath package (from 
the AMS Euler fonts). If you don't have the AMS font set, you will have to use the 
normal math italic symbols and the typesetter will substitute the corresponding 
upright characters. 

You will make this much easier if you can use the macros \upi and \umu in your 
text to indicate the need for the upright characters, together with the following 
definitions in the preamble, where they can easily be redefined by the typesetter. 

\providecoinincind\upi{\pi} 
\providecoinmajid\umu{\mu)- 

Notice the use of \providecommand, this stops your definitions overwriting the 
ones used by the typesetter. 

Sloping uppercase Greek characters can be obtained by using the \varGamma etc. 
macros (which are built into the class file). The few required upright Greek capital 
letters (e.g. capital delta for 'difference') are given by the relevant X^T^jX macro, 
e.g. \Delta. 
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4-3.4- Sans serif symbols The \textsf and \mathsf commands change the type- 
face to sans serif, giving upright characters. Occasionally, bold-sloping sans serif is 
needed. You should use the following supplied macros to obtain these fonts. 

\textsf {text} — !■ text \mathsf {math} — > math 

\textsf i{text} — > text \mathsf i{math} — > math 

\textsfb{text} — > text \mathsfb{inath} — > math 

\textsfbi{text} — > text \inathsfbi{inath} — > math 

You can use them like this: 

\newcoinmcind\ssC{\mathsf {C}} 7. for sans serif C 
\newcoinmand\sf sP{\mathsf i{P}} % for sans serif slanted P 
\newcoinmand\sfbsX{\inatlisfbi{X}} % for sans serif bold slanted X 

Note that the bold-slanted font \testsf bi and \mathsf bi use the slanted sans serif 
font cms si because there is no bold-slantcd maths sans serif font in Computer 
Modern! If you use the supplied sans-serif text and math commands the typesetter 
will be able to substitute the fonts automatically. 

4-3.5. Bold math italic symbols If you require bold math italic symbols/letters, 
provides several ways of getting them. Firstly, X^T^jK's \boldmath switch 
which can be used like this: 

$ \mbox{\boldmath $P$} = \mathsf {J}\mathcal{E} T $ 

As you can see it takes quite alot of typing to achieve a bold math italic P. Another 
problem is that you can't use \boldmath in math mode, thus the \mbox forces the 
resulting text into text style (which you may not want). 

You can cut down on the typing by defining the following macro (which is in the 
preamble of this guide): 

\providecommand\boldsymbol [1] {\mbox{\boldmath $#1$}} 

Then for the above example you can define: 

\newcommand\bldP{\boldsymbol{P}} 

the achieve the same result. This still doesn't remove the default to text style 
problem. 

If you have the amsbsy package on your system, then you can remove this lim- 
itation as well by placing a \usepackage{amsbsy} after the \dociimentclass line 

in your document, and then use \boldsymbol for bold math italic symbols (don't 
define \boldsymbol yourself when using cimsbsy). 

4-3-6- Script characters Script characters should be typeset using 2e 's \mathcal 
command. This produces the Computer Modern symbols such as £ and in your 
hard copy but the typesetter will substitute the more florid script characters nor- 
mally seen in the journal. 

4-4- Skewing of accents 

Accents such as hats, over bars and dots are normally centred over letters, but when 

these are italic or sloping greek the accent may need to be moved to the right so that 

it is centred over the top of the sloped letter. For example, \newcommand\hatp{\skew3\hat{p}} 

will produce p. 
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4-5. Units of measure 

Numbers and their units of measure should be typeset with fixed spaces that will 
not break over two lines. This is easily done with user-defined macros. For example, 
52\dynpercm typesets as 52 dynes cm~^, providing the following macro definition 

has been included in the preamble. 

\newcoininEind\dynpercm{\nobreak\mbox{$\ ; $dynes\ , cm$~{-l}$}3- 

5. User-defined macros 

If you define your own macros you must ensure that their names do not conflict 
with any existing macros in plain T^TjX or I^T[t;X (or AMSI^T[t;X if you are using 
this). You should also place them in the preamble to your input file, between the 
\dociimentclass and \begin-C document} commands. 

Apart from scanning the indexes of the relevant manuals, you can check whether 
a macro name is already used in plain TjtjX or I^T^t;X by using the command 
show. For instance, run l^T^X interactively and type \show\<macro_name> at the 
T^pC prompt. (Alternatively, insert the \show command into the preamble of an 
input file and T^TjX it.) 

* \show\Re 

produces the following response from TJn^X: 

> \Re=\mathchar"23C. 
<*> \show\Re 

By contrast \Real is not part of plain T^jX or I^T^jK and \show\Real generates: 

> \Real=undef ined 
<*> \show\Real 

confirming that this name can be assigned to a user-defined macro. 

Alternatively, you can just use Xnewcommand, which will check for the existence 
of the macro you are trying to define. If the macro exists will respond with: 

! LaTeX Error: Command ... already defined. 

In this case you should choose another name, and try again. 

Such macros must be in a place where they can easily be found and modified 
by the journal's editors or typesetter. They must be gathered together in the 
preamble of your input file, or in a separate macros.tex file with the command 
\input -[macros} in the preamble. Macro definitions must not be scattered about 
your document where they are likely to be completely overlooked by the typesetter. 

The same applies to font definitions that are based on Computer Modern fonts. 
These must be changed by the typesetter to use the journal's TJ^X typefaces. In 
this case, you should draw attention to these font definitions on the hard copy that 
you submit for publication and by placing a comment in your input file just before 
the relevant definitions, for example % replace font! 

6. Some guidelines for using standard facilities 

The following notes may help you achieve the best effects with the standard M^iX 
facilities that remain in the JPP style. 
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Table 1. An example table 
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— 2 
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6.1. Sections 

provides five levels of section headings, only four of which are defined in the 
JPP class file: 

Heading A - \section{. . .} 
Heading B - \subsection{. . .} 
Heading C \subsubsection-[ . . .} 
Heading D \paragraph{. . .} 
There is no subparagraph in the JPP style. 

To obtain non-bold in a bold heading use the usual WTf^2s commands for 
changing typeface; for example \section{Fluctuations in Ca\textsc{ii». 

6.2. Tables 

The table environment is implemented as described in the I^T^riX manual to provide 
consecutively numbered floating inserts for tables. 

JPP class will cope with most table positioning problems and you should not 
normally use the optional positional qualifiers t, b, h on the table environment, 
as this would override these decisions. Table captions should appear at the bot- 
tom of the table; therefore you should place the \caption command before the 
\begin-[tabular>. 

The JPP table environment will insert rules above and below the table as re- 
quired by the JPP style, so you should not attempt to insert these yourself. The 
only time you need to intervene with the rules is when two or more tables fall one 
above another. When this happens, you will get two rules together and too much 
space between the tables; the solution is to add \norule before the end of upper 
table and \f ollowon between tables as shown here: 

\end{center} 

\caption{Dimensionless paraineters}\label-Cdimensionless_p} 

\norule 
\end{table} 

\f ollowon 
\begin{table} 

\begin{center} 

The JPP style dictates that vertical rules should never be used within the body of 
the table, and horizontal rules should be used only to span columns with the same 
headings. Extra space can be inserted to distinguish groups of rows or columns. 

As an example, table 1 is produced using the following commands: 

\begin{table} 

\caption-[An example table} \label-[sample-table} 
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\begin{center} 

\begin-(tabular}-{cccc} 
\hline 

{Figure} & {$hA$} & ■C$hB$} & {$hC$}\\ 
\hline 

2 & $\exp\;(\upi \mathrm{i} x)$ 

& $\exp\; (\iipi \mathrm{i} y)$ & $0$\\ 

3 & $-1$ & $\exp\;(\upi \mathrm{i} x)$ & $1$\\ 

4 & $-4+3\matlirm{i}$ & $-4+3\mathrm{i}$ & 1.6\\ 

5 & $-2$ & $-2$ & $1.2 \mathrm{i}$ 
\end-[tabular}- 

\end{center]- 
\end{table} 

The tabular environment has been modified for the JPP style in the foUowing 
ways: 

(a) Additional vertical space is inserted above and below a horizontal rule produced 
by Xhline 

(6) Tables are centred, and span the full width of the page; that is, they are similar 
to the tables that would be produced by \begin{tabular*}{\textwidth}. 

(c) As with normal I^li^X, if a footnote needs to be inserted into the body of the 
table, the tabular environment should be enclosed in a minipage environment 
of width \textwidth. e.g. 

\begin{minipage}{\textwidth} 
\begin-Ctabular}-[cccc}- 
Xhline 

{Figure} & {$hA$} & {$hB$}\f ootnote{A table must be 
inside a \verb"ininipage" environment if it includes 
table footnotes.} 
& {$hC$}\\ 

\end{tabular} 
\end{minipage} 

Commands to redefine quantities such as \arraystretch should be omitted. If 
the original tabular facilities are needed, there is a new environment 'oldtabular', 
which has none of the reformatting; it should be used in exactly the same way. 

6.3. Illustrations (or figures) 

Artwork for all figures must be supplied separately, to be processed by the Printer, 
and so I^TJt;X's picture environment, for example, cannot be used. An approximate 
amount of space should be left, using the \vspace command. 

The JPP style will cope with most figure positioning problems and you should not 
normally use the optional positional qualifiers t, b, h on the figure environment, 
as this would override these decisions. Figure captions should be below the figure 
itself, therefore the \caption command should appear after the space left for the 
illustration within the figure environment. For example, figure 1 is produced using 
the following commands: 
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Figure 1. An example figure with space for artwork. 

\begin{f igure} 
\vspace{3cm} 

\caption{An example figure with space for artwork. } 
\label{sample-f igure} 
\end{f igure} 

If a figure caption is too long to fit on the same page as its illustration, the caption 

may be typeset as 'Figure X. For caption sec facing page.', and the longer caption 
typeset at the bottom of the facing page. Authors should not concern themselves 
unduly with such details, and may leave pages long. 

6.4- Acknowledgments 

Acknowledgments should appear at the close of your paper, just before the list 
of references and any appendices. Use the acknowledgments environment, which 
simply inserts an 'Acknowledgments' subsection heading. 

There is also an acknowledgment (or acknowledegment) environment for type- 
setting a single acknowledegment. 

6. 5. Appendices 

You should use the standard WI^^ \appendix command to place any Appendices, 
normally, just before the references. This numbers appendices as A, B etc. , equa- 
tions as (Al), (Bl) etc. , and figures and tables as Al, Bl etc. 

If you have only one Appendix, you should use the \oneappendix command 
(instead of \appendix). This command does the same as \appendix, except it 
allows jpp.cls to typeset a single Appendix correctly. 

6. 6. References 

As with standard there are two ways of producing a list of references; either 

by compiling a list (using a thebibliography or theref erences environment), or 

by using BibTjrjX with a suitable bibliographic database. 

References in the text should follow either the Harvard (name/date) or Vancouver 
(numbered) system. In the Harvard system, each reference should be cited in the 
text as 'author(s) (date)' or '(author(s) date)', and the reference list at the end 
of the paper should be in alphabetical order. In the Vancouver system, references 
should be indicated in the text by numbers in square brackets (e.g. [1], [1,2], [1-4], 
etc.), and should be listed at the end of the paper in numerical order of citation. 
For details of the style of entries within the reference list, please see any recent issue 
of the journal. 

6.6.1. References in the text References in the text are given by author and date 
like Dennis (1985) or numbered system [12]. Each entry has a key, which is assigned 
by the author and used to refer to that entry in the text. 
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6.6.2. Harvard (name/date) The following listing shows some references prepared 
in author (date) style of the journal; the code produces the references at the end of 
this guide. 

\begin{thebibliography}{} 

\bibitein[Abramowitz and Stegun (1965)] {AS65} 
{Abramowitz, M. and Stegun, I.\,A.> 1965 
{\ein Handbook of Mathematical Functions}. Dover. 

Xbibitem [Dennis (1985)] {Den85} 
{Dennis, S.\,C.\,R.} 1985 

Compact explicit finite difference approximations to the 
Navier — Stokes equation. In {\em Ninth Intl Conf . on 
Numerical Methods in Fluid DynaimicsX/} (ed. Soubbaramayer 
and J.\,P. Boujot) . Lecture Notes in Physics, vol. 218, 
pp. 23 — 51. Springer. 
\bibitem [Jones (1976)] {Jonl976} 
{Jones, 0.\,C.} 1976 

An improvement in the calculation of turbulent friction in 
rectangular ducts. {\em Trans. ASME\/} J: 
{\em J. \, Fluids Engng\/> {\bf 98}, 173—181. 
\bibitem [Saf f man ( 1990) ] {Saf 1990} 
{Saffman, P.\,G.} 1990 

A model vortex reconnection. {\em J. \, Fluid Mech.} 

{\bf 212>, 395—402. 
\bibitem [Saffman and Schatzman (1982)] {SS1982} 

{Saffman, P.\,G. and Schatzman, J.\,C.} 1982 

Stability of a vortex street of finite vortices. 

{\em J. \, Fluid Mech.} {\bf 117}, 171—185. 
\bibitem [Saffman and Yuen (1980)] {SY1980} 

{Saffman, P.\,G. and Yuen, H.\,C.} 1980 

A new type of three-dimensional deep-water wave of permanent 
form. {\em J. \, Fluid Mech.} {\bf 101}, 797—808. 

\bibitem[Shaqfeh and Koch (1990)] {SK1990} 
{Shaqfeh, E.\,S.\,G. and Koch, D.\,L.} 1990 
Orientational dispersion of fibres in extensional flow. 
{\em Phys. Fluids\/} A {\bf 2}, 1077—1081. 

\bibitem [Wi jngaarden (1968) ] {Wi j 1968} 
{Wijngaarden, L. van} 1968 

On the oscillations near and at resonance in open pipes. 
{\em J.\,Engng Maths\/} {\bf 2}, 225—240 
\bibitem [Williams (1964)] {Wil64} 
{Williams, J.\,A.} 1964 

A nonlinear problem in surface water waves. PhD thesis. 
University of California, Berkeley. 
\end{thebibliography} 

Each entry takes the form 



\bibitem [Author (s) (Date)] {tag} 
Bibliography entry 
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where Author (s) should be the author names as they are cited in the text, Date is 
the date to be cited in the text, and tag is the tag that is to be used as an argument 
for the \cite{} command. Bibliography entry should be the material that is to 
appear in the bibliography, suitably formatted. 

References 

Abramowitz, M. and Stogun, I. A. 1965 Handbook of Mathematical Functions. Dover. 
Dennis, S. C. R. 1985 Compact explicit finite difference approximations to the Navier- 

Stokes equation. In Ninth Intl Conf. on Numerical Methods in Fluid Dynamics (ed. 

Soubbaramayer and J. P. Boujot). Lecture Notes in Physics, vol. 218, pp. 23-51. 

Springer. 

Jones, O. C. 1976 An improvement in the calculation of turbulent friction in rectangular 

ducts. Trans. ASME J: J, Fluids Engng 98, 173-181. 
Saffman, P. G. 1990 A model vortex reconnection. J. Fluid Mech. 212, 395-402. 
Saffman, P. G. and Scliatzman, J. C. 1982 Stability of a vortex street of finite vortices. 

J. Fluid Mech. 117, 171-185. 
Saffman, P. G. and Yuen, H. C. 1980 A new type of three-dimensional deep-water wave of 

permanent form. J. Fluid Mech. 101, 797-808. 
Shaqfeh, E. S. G. and Koch, D. L. 1990 Orientational dispersion of fibres in extensional 

flow. Phys. Fluids A 2, 1077-1081. 
Wijngaarden, L. van 1968 On the oscillations near and at resonance in open pipes. J. Engng 

Maths 2, 225-240 

Williams, J. A. 1964 A nonlinear problem in surface water waves. PhD thesis. University 
of California, Berkeley. 

6.6.3. Vancouver (numbered) system The following listing shows some references 
prepared in Vancouver (numbered) system style of the journal; the code produces 
the references at the end of this guide. 

\begin{theref erences}-[9} 
\bibitem{Jon76} 
{Jones, 0.\,C.}, 

An improvement in the calculation of turbulent friction in 
rectangular ducts. {\em Trans. ASME\/} J: 
{\em J. \, Fluids Engng\/> i\M 98}, 173—181 (1976). 
\bibitem{Saf90} 
{Saffman, P.\,G.>, 

A model vortex reconnection. {\em J. \, Fluid Mech.} 

{\bf 212}, 395—402 (1990). 
\bibitem{SS82} 

{Saffman, P.\,G. and Schatzman, J.\,C.}, 

Stability of a vortex street of finite vortices. 

{\em J. \, Fluid Mech.} {\bf 117}, 171—185 (1982). 
\bibitem{SY80} 

{Saffman, P.\,G. and Yuen, H.\,C.}, 

A new type of three-dimensional deep-water wave of permanent 
form. {\em J. \, Fluid Mech.} {\bf 101}, 797—808 (1980). 
\bibitem{SK90} 

{Shaqfeh, E.\,S.\,G. and Koch, D.\,L.} 

Orientational dispersion of fibres in extensional flow. 
{\em Phys. Fluids\/} A {\bf 2}, 1077—1081 (1990). 
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\bibitem-CWij68} 

{Wijngaarden, L. van}, 

On the oscillations near and at resonance in open pipes. 
{\em JA.Engng Maths\/} {\bf 2}, 225—240 (1968). 
\end{theref erences} 

Each entry takes the form 

Xbibitem-Ctag} 

Bibliography entry 

where tag is the tag that is to be used as an argument for the \cite{} command. 
Bibliography entry should be the material that is to appear in the bibliography, 
suitably formatted. 

References 

[1] Jones, O.C., An improvement in the calculation of turbulent friction in rect- 
angular ducts. Trans. ASME J: J. Fluids Engng 98, 173 181 (1976). 

[2] Saffman, P. G., A model vortex reconnection. J. Fluid Mech. 212, 395-402 
(1990). 

[3] Saffman, P. G. and Schatzman, J. C., Stability of a vortex street of finite vor- 
tices. J. Fluid Mech. 117, 171-185 (1982). 

[4] Saffman, P. G. and Yuen, H. C., A new type of three-dimensional deep-water 
wave of permanent form. J. Fluid Mech. 101, 797-808 (1980). 

[5] Shaqfeh, E. S. G. and Koch, D. L. Orientational dispersion of fibres in exten- 
sional flow. Phys. Fluids A 2, 1077-1081 (1990). 

[6] Wijngaarden, L. van, On the oscillations near and at resonance in open pipes. 
J. Engng Maths 2, 225-240 (1968). 
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Appendix A. Special commands in jpp.cls 

The following is a summary of the new commands, optional arguments and envi- 
ronments which have been added to the standard IM^jX user-interface in creating 
the JPP class file. 



New commands 

\af filiation 
\af f ilskip 

\author 
\f ollowon 
\ls 

\newtheorem 



\newdef inition 



\norule 
\ns 

\oneappendix 

New optional arguments 
[<short title>] 

[<short author>] 



[<widest label>] 

New environments 

acknowl edg ( e ) ment ( s ) 

proof 

proof* 

remark 



subeqnarray 
subequation 
tabular 



use after \author to typeset an affiliation, 
use immediately after \\ in \af filiation to give 
correct vertical spacing between separate affiliations, 
do not use \\ in \author to start an affiliation, 
to remove space between adjacent tables, 
to letter space the author's name, 
this is enhanced so that you can produce unnum- 
bered versions of the environments by using the * 
form. e.g. \begin{theorem*}. 

this environment works like the Xnewtheorem com- 
mand, except it produces environments which are 
typeset in roman instead of italic. Again unnumbered 
'definitions' can be typeset using the * forms, 
to remove the rule below a table, 
to insert space between an author's names, 
does the same as \appendix, except it allows 
jpp. els to typeset a single Appendix correctly. 



in the \title command: to define a right running 
headline that is different from the article title, 
in the \author command: to define a left running 
headline that is different from the authors' names as 
typeset at the article opening. 

in \begin{ enumerate}: to ensure the correct align- 
ment of numbered lists. 



to typeset acknowledgements, 
to typeset mathematical proofs, 
to typeset mathematical proofs without the termi- 
nating proofbox. 

this environment works like the theorem environ- 
ment; it typesets an italic heading and roman text 
to contrast with theorem's small caps heading and 
italic text. 

enables equation numbers in an array to be num- 
bered as (6.1a), (6.2b), etc. 

enables consecutive equations to be numbered (6.1a), 
(6.1b), etc. 

has been modified to insert additional space above 
and below an \hrule and the table caption and body 
is centred with rules full out across the text measure. 
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Appendix B. Notes for editors 

This appendix contains additional information which may be useful to those who are 
involved with the final production stages of an article. Authors, who are generally 
not typesetting the final pages in the journal's typeface (Monotype Modern), do 
not need this information. 

B.0.1. Catchline commands To be placed in the preamble: 

• \date{> 

• \pagerange{} 

• \part{} 

• \pubyear{} 

• \volume{} 

• \doi{> 

B.0.2. Footnotes If a footnote falls at the bottom of a page, it is possible for the 
footnote to appear on the following page (a feature of TJ^X). Check for this. 

B.0.3. Rules between tables If two or more tables fall one above another, add the 
commands \norule and \f ollowon as described in this guide. 

B.0.4- Font substitution Check for use of AMS fonts, bold slanted sans serif, and 

bold math italic and alter preamble definitions to use the appropriate AMS / CUP /Monotype 

fonts for phototypesetter output. 

B.0.5. Font sizes The JPP class file defines all the standard font sizes. For 

example: 

• tiny — This is tiny text. 

• SCriptsize — This is sciiptsizc text. 

• \indexsize — This is indexsize text. 

• \f ootnotesize - This is footnotesize text. 

• \ small - This is small text. 

• \qsmall - This is qsmall text (quotations). 

• \normalsize - This is normalsize text (default). 

• Marge - This is large text. 

• \Large - This is Large text. 

• \LARGE- This is LARGE text. 

All these sizes are summarized in Table B.l. 
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Table B.l. Type sizes for M^jK size-changing commands 



Size 


Size/Baseline 


Usage 


\tiny 


5pt/6pt 




\sc!rip"t siz6 


/ pi / opi 




\indexsize 


8pt/9pt 


catchline. 


\f ootnotesize 


9pt/10pt 


index, footnotes. 


\ small 


9pt/10pt 


'AND' in authors' names, re- 






ceived date, affihation, fig- 






ure and table captions. 


\qsmall 


9.75pt/10.75pt 


quote, quotations. 


Xnormalsize 


10pt/12pt 


main text size, abstract, 




B and C heads. 


Marge 


llpt/13pt 


A head, part no. 


\Large 


14pt/18pt 


part title. 


\LARGE 


17pt/19pt 


article title. 


\liuge 


20pt/25pt 




\Huge 


25pt/30pt 
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An interesting and seminal work on 
various phenomena in fluid mechanics 

ALAN N. OTHERH, 
H.-C. SMITHi and J.Q. PUBLIC^ 

^Department of Chemical Engineering, University of America, Somewhere, IN 12345, 

USA 

^Department of Aerospace and Mechanical Engineering, University of Camford, 
Academic Street, Camford, CF3 5QL, UK 

(Received 1 July 1994 and in revised form 30 October 1998) 

Abstract. Using Stokes flow between eccentric counter-rotating cylinders as a 
prototype for bounded nearly parallel lubrication flow, we investigate the effect of 
a slender recirculation region within the flow field on cross-stream heat or mass 
transport in the important limit of high Peclet number Pe where the enhancement 
over pure conduction heat transfer without recirculation is most pronounced. The 
steady enhancement is estimated with a matched asymptotic expansion to resolve 
the diffusive boundary layers at the separatrices which bound the recirculation 
region. The enhancement over pure conduction is shown to vary as at infinite Pe, 
where is the characteristic width of the recirculation region. The enhancement 
decays from this asymptote as Pe ^ . 



1. Introduction 

The use of integral equations to solve exterior problems in linear acoustics, i.e. to 
solve the Helmholtz equation (V^ -f = outside a surface S given that 

satisfies certain boundary conditions on S", is very common. A good description is 
provided by 9. Integral equations have also been used to solve the two-dimensional 
Helmholtz equation that arises in water-wave problems where there is a constant 
depth variation. The problem of wave oscillations in arbitrarily shaped harbours 
using such techniques has been examined (see for example Hwang & Tuck 1970; 
Lee 1971). 

In a recent paper 8 have shown how radiation and scattering problems for vertical 
circular cylinders placed on the centreline of a channel of finite water depth can be 
solved efficiently using the multipole method devised originally by Ursell (1950). 
This method was also used by 2 to prove the existence of trapped modes in the 
vicinity of such a cylinder at a discrete wavenumber k < 7r/2(i where 2d is the 
channel width. 

Many water- wave/body interaction problems in which the body is a vertical 
cylinder with constant cross-section can be simplified by factoring out the depth 
dependence. Thus if the boundary conditions are homogeneous we can write the 
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velocity potential (j){x,y,z,t) = Re{(j){x,y) cosh k{z + /i)e~"^*}, where the {x,y)- 
plane corresponds to the undisturbed free surface and z is measured vertically 
upwards with z = —h the bottom of the channel. 

Subsequently Callan et al. (1991) proved the existence of, and computed the 
wavenumbers for, the circular cross-section case. It should be noted however that 
experimental evidence for acoustic resonances in the case of the circular cylinder is 
given by Bearman & Graham (1980, pp. 231-232). 

6 provided a theory for determining the trapped-mode frequencies for the thin 
plate, based on a modification of the Wiener-Hopf technique. 

The use of channel Green's functions, developed in § 2, allows the far-field be- 
haviour to be computed in an extremely simple manner, whilst the integral equation 
constructed in § 3 enables the trapped modes to be computed in § 4 and the scatter- 
ing of an incident plane wave to be solved in § 5. Appendix A contains comparisons 
with experiments. 



2. Green's functions 

2.1. Construction of equations 

We are concerned with problems for which the solution, 0, is either symmetric or 
antisymmetric about the centreline of the waveguide, y = 0. The first step is the 
construction of a symmetric and an antisymmetric Green's function, Gs{P, Q) and 
Ga{P,Q)- Thus we require 

(Vl + fc2)G« = (V| + k^)Ga = (2.1) 

in the fluid, where Vp is a gradient operator, except when P = Q, 

Gs,Ga~l/(27r)lnr asr=|P-Q|^0, (2.2) 

dG, dG, 



dy dy 
dG, 



on y = d, (2.3) 



= ony = 0, (2.4) 
dy 

Ga on y = 0, (2.5) 

and we require Gg and Ga to behave like outgoing waves as |a;| — > oo. 
One way of constructing Gs or Ga is to replace (2.1) and (2.2) by 

(V^ + fc2)G, = (V^ + fc2)Ga = -5{x - ^)S{y - rj) (2.6) 

and to assume initially that k has a positive imaginary part. 

2.2. Further developments 

Using results from 8 we see that this has the integral representation 

1 poo 

/ [e-'^'^l^-"' -t- e-'^^^^d-y-r,)! _ Q^y^ ^ ^ (2.7) 

27r Jo 

where 

'^'^ = \ it^-l)K t>l. 
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In order to satisfy (2.4) wc add to this the function 
1 f°° ^^^^ cosh fc7(d — y) 

which satisfies (2.1), (2.3) and obtain 

B{t) = 2e-'=^'* cosh kj{d - ??). (2.8) 

Thus the function 

1 f°° Q-^ld, 

G = -^i{Ho{kr)+Ho{kn))-- j- . ^, ^ cosh fc7(rf-y) cosh fc7(rf- 77) (2.9) 

TT Jq J sinh k^d 

satisfies (2.1)-(2.4). By writing this function as a single integral which is even in 7, 
it follows that G is real. 



3. The trapped-mode problem 

The unit normal from D to dD is = {-y'{0),x'{9))/w{9). Now G„ = jYo{kr) + 
Ga where r = {[x{6) - x{ip)]'^ + [y{6) — J/(V')]^}^ and Ga is regular as kr — »■ 0. In 
order to evaluate dGa{0, 9)/dnq we note that 

as kr —>■ 0. Expanding x{^jJ) and y{ip) about the point ip = 9 then shows that 

^ -[p'{0)p"{e)-p'^{e)-2p'^{e)] askr^O. (3.1) 



477^3 (6»)' 
5.-?. Computation 

For computational purposes we discretize (3.1) by dividing the interval (0,7r) into 



M segments. Thus we write 

M 



^'^W = ^E'^(^i)a^^«(^'^^X^^)' 0<V'<7r, (3.2) 
where 0^- = (j — ^)it/M. Collocating at ip = 9i and writing = (/'(^i) etc. gives 

M 

^•^^ = ^ E '^i^S^'j' i = 1, • • • , M, (3.3) 



where 

1 aG„(e,, e,)/an, + [p^p^ -Pi- 2pf]/A-Kwl i = j. ^ - ' 

For a trapped mode, therefore, we require the determinant of the M x M matrix 
whose elements are 

27r 

to be zero. Table 1 shows a comparison of results obtained from this method using 
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Table 1. Values of kd at which trapped modes occur when p{0) = a 



a/d 


M = 4 


M = 8 


Callan et al. 


0.1 


1.56905 


1.56905 


1.56904 


0.3 


1.50484 


1.50484 


1.50484 


0.5 


1.39128 


1.39131 


1.39131 


0.7 


1.32281 


1.32287 


1.32288 


0.9 


1.34479 


1.35160 


1.35185 



Figure 1. Trapped-mode wavenumbers, kd, plotted against a/d for three ellipses: 
, b/a = 0.5; , b/a = 1; , b/a = 1.5. 

two different truncation parameters with accurate values obtained using the method 
of Callan et al. (1991). 

An example of the results that are obtained from our method is given in figure 1. 

Figure 2(a,b) shows shaded contour plots of (j) for these modes, normalized so that 
the maximum value of (j) on the body is 1. 

3.2. Basic properties 
Let 

pi= lim p{x,C}, pu = lim p{x,C) (3.5a,b) 

be the fluid densities immediately below and above the cat's-eyes. Finally let po 
and A^o be the constant values of the density and the vorticity inside the cat's-eyes, 
so that 

{p{x,0,(l>aix,0) = {po,No) ioT Zlix) < C < Z^{x). (3.6) 
The Reynolds number Re is defined by Ut-H/v {y is the kinematic viscosity), the 
length given in wall units is denoted by ( )+, and the Prandtl number Pr is set 
equal to 0.7. In (2.1) and (2.2), and are 

Tij = [uiUj - UiUj) + {uiU^'^^ + uf^^Uj) + uf'^^uf'^^, (3.7a) 

= (^ - Uje) + {UjeSGS + u^GSQ-^ + uSGS0SGS_ (3 Jb) 
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Figure 2. Sliaded contour plots of the potential (f> for the two trapped modes that exist 
for an ellipse with a/d = 1.5, b/d = 0.75. (a) Symmetric about x = 0, kd = 0.96; (6) 
antisymmetric about x = 0, kd = 1.398. 

3.2.1. Calculation of the terms The first terms in the right-hand side of (3.5a) 
and (3.5b) are the Leonard terms exphcitly calculated by applying the Gaussian 
filter in the x- and ^-directions in the Fourier space. 

The interface boundary conditions given by (2.1) and (2.2), which relate the dis- 
placement and stress state of the wall at the mean interface to the disturbance 
quantities of the flow, can also be reformulated in terms of the transformed quanti- 
ties. The transformed boundary conditions are summarized below in a matrix form 
that is convenient for the subsequent development of the theory: 



Qc = 



o w 



iV(a*+u;-iO 



i/3 



Rx V'' 



-(a*a;)-i 






(3i?,-i+c*(ia*)-i) 








-(ia*i?5)-i 















(3.8) 

is termed the displacement-stress vector and Qc the flow-wall coupling matrix. 
Subscript w in (3.8) denotes evaluation of the terms at the mean interface. It is 
noted that = for the Blasius mean flow. 

3.2.2. Wave propagation in anisotropic compliant layers From (2.1), the funda- 
mental wave solutions to (3.1) and (3.2) for a uniformly thick homogeneous layer 
in the transformed variables has the form of 



Tj* = T7* exp[i(a*x* - cot)], 



(3.9) 



where r)* ~ 6cxp(i7a:;| 
and (1.4) yields the following detcrminantal equation: 



For a non-trivial wave, the substitution of (1.7) into (1.3) 



(3.10) 
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where the wavenumbers {k{, k^, k^) = (a*, 0, 7) and 5ps is the Kronecker delta. 



4. Torus translating along an axis of symmetry 

Consider a torus with axes a,b (sec figure 2), moving along the z-axis. Symme- 
try considerations imply the following form for the stress function, given in body 
coordinates: 

f{0, V) = (5(V) cos e, (V) sin e, /(V)). (4.1) 

Because of symmetry, one can integrate analytically in the ^-direction obtaining 
a pair of equations for the coefficients /, g in (4.1), 

''^'^ (sin -01 — sin-0)(a + bcosip)i 



/(V-i) 



3b 



fi'(V'i) 



7r(2(a + 6cosV'i)) 
3 



r'' (sin7/> 

i))iJo (1- 



(4.2) 



'ir{2{a + bcosipi))2 Jo 
sin ijji — sin tp 



cos(V'i — V') 



^'^ ^a + b cos ip 



(1 - cos(V'i - V))(2 + a)2 

/(V')[(cos^i - cos V - bPi)S + /3iP] 



2 + a 

(sinV^i -sinV>)^ 2 , o 

2 + a- -— - 6 7 5 

1 — cos(V' — V'l) 



+ 



(b^ cos V'17 - ^aj F^TT, ^) - (2 + a) costpiE{^'K, 6) 

^2(1 -COS(V'-V'I)) 



(a + fa cos '(/^) (a + facos'i/'i) ' 



/3-/?(^,Vi) = 



dV, (4.3) 

1 — C0S('(/) — t/ji) 

a + b cos ip 

(4.4) 



5. Conclusions 

We have shown how integral equations can be used to solve a particular class of 

problems concerning obstacles in waveguides, namely the Neumann problem for 
bodies symmetric about the centreline of a channel, and two such problems were 
considered in detail. 
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Appendix A. Boundary conditions 

It is convenient for numerical purposes to change the independent variable in (4.1) 

1 1 

to z = y/v^ and to introduce the dependent variable H{z) = (/ — y)/v^. Equation 
(4.1) then becomes 

(1 -P)iH + z)H" - (2 + H')H' = H'". (A 1) 
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Boundary conditions to (4.3) follow from (4.2) and the definition of H: 
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if'(0) = -l + eiC7„ + eC4; 




H'{oo) = 0. 



> 



(A 2) 




Appendix B 

A simple sufficient condition for the method of separation of variables to hold for 
the convection problem is derived. This criterion is then shown to be satisfied for 
the ansatz described by (3.27), thus justifying the approach used in §3. The basic 
ingredient of our argument is contained in the following estimate for a Rayleigh 
Ritz ratio: 

Lemma 1. Let f{z) be a trial function defined on [0, 1]. Let Ai denote the ground- 
state eigenvalue for — d^g/dz^ = Ag, where g must satisfy ±dg/dz + ag = at 
z = 0,1 for some non-negative constant a. Then for any f that is not identically 
zero we have 



Before proving it, we note that the first inequality is the standard variational 
characterization for the eigenvalue Ai . 

Corollary 1. Any non-zero trial function f which satisfies the boundary con- 
dition /(O) = /(I) = always satisfies 
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Abstract. 

We perform a fully self-consistent 3-D numerical simulation for a compressible, 
dissipative magneto-plasma driven by large-scale perturbations, that contain a fairly 
broader spectrum of characteristic modes, ranging from largest scales to interme- 
diate scales and down to the smallest scales, where the energy of the system are 
dissipated by coUisional (Ohmic) and viscous dissipations. Additionally, our sim- 
ulation includes nonlinear interactions amongst a wide range of fluctuations that 
are initialized with random spectral amplitudes, leading to the cascade of spectral 
energy in the inertial range spectrum, and takes into account large scale as well as 
small scale perturbation that may have been induced by the background plasma 
fluctuations, also the non adiabatic exchange of energy leading to the migration of 
energy from the energy containing modes or randomly injected energy driven by 
perturbations and further dissipated by the smaller scales. Besides demonstrating 
the comparative decays of total energy and dissipation rate of energy, our results 
show the existence of a perpendicular component of current, thus clearly confirming 
that the self-organized state is non-force free. 



1. Introduction and Motivation: 

The phenomena of self-organization, in which a continuous system naturally evolves 
towards a state exhibiting some form of order on large scales, is deeply rooted in 
nature. The ordered states are remarkably robust; their detailed structure remains 
relatively invariant across experimental realization; these preferred states are inde- 
pendent of the way the system is prepared (Hasegawa, 1985; Ortolani and Schnack, 
1993) 

Several broad classes of physical processes , such as first and second order phase 
transitions, crystallization, structure formations in space and astrophysics and cos- 
mology can be described as examples of self-organization. In chemistry, examples of 
self-organizations include reaction-diffusion system such as Belousov-Zhabotinsky 
reactions, self- assembled monolayers, molecular self-assembly, etc. Biological sys- 
tems are swarmed with examples of self-organizations; besides the example of pat- 
tern formation and morphogenesis, the origin of life itself from self-organizing chem- 
ical systems is the supreme example self-organization. Most of the system exhibiting 
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self-organization in nature are intrinsically nonlinear and often are not isolated - 
they are driven by some external input. Self-organization occurs in open systems 
which are far from thermal equilibrium. Nicolis and Prigogine (1977) had envisaged 
new types of self-organized states for driven dissipative systems and called these 
states as "Dissipative Structures". Such structures provide striking examples of 
non-equilibrium as a source of order. During the past decades, non-linear dynamics 
and non- equilibrium thermodynamics have developed along with plasma physics, 
providing qualitatively new approaches to complex problems. It is thus widely ac- 
cepted nowadays, for instance, that non-equilibrium may be a source of order in 
dissipative systems, allowing the emergence of self-organization Self-organizations 
have been observed in varieties of laboratory magnetically confined systems. Most 
plasma systems are dissipative. externally driven and far from equilibrium, thus 
sharing some essential features with other complex non-linear systems that show 
spatial and temporal coherence. Such a system is described by a set of non-linear 
partial differential equations. During its relaxation towards a self-organized state, it 
takes advantage of instabilities that lead it, for instance, to a preferred state, under 
certain constraints that prevent it from falling into a trivial unconfined state. Such 
constraints appear as quadratic or higher order quantities, which are conserved in 
the absence of dissipation. The relevant feature is that, in the presence of dissipa- 
tion one of these quantities decays faster than the others. Generally, it remains to 
be determined which is the relevant variational principle underlying the relaxation 
mechanisms described above, is appropriate to characterize the self-organized state. 
In the space and solar plasma flows, a number of commonly observed processes in- 
cluding flares, prominences, filaments and/or (magnetic) loop like structures that 
are generated often during the active solar period can be modelled by relaxation 
of plasma through self organization (Bhattacharya et al., 2007). The kinematic as 
well as dynamics of these entities are far more complex than ever thought and we 
still lack an in depth insight into their evolutionary characteristics despite the avail- 
ability of an enormous data bases from various spacecraft missions. A part of the 
problem lies with the lack of a fully self-consistent description (or physical model) 
of evolution of these structures in realistic environments. For instance, one of the 
mechanisms of the generation of magnetic field loops is often attributed to the Tay- 
lor relaxation (Taylor, 1974) process where magnetic stresses overcome the pressure 
stresses which thereby balance all the MHD forces. Such state is often character- 
ized by a low plasma beta (where plasma-beta is a ratio of pressure and magnetic 
energy). Under no external forces, the pressure gradients flatten out to nullify the 
magnetic pinch (J x B) force. This state is called a "Force Free State". The force 
free state is one of the analytic descriptions that describes magnetic loop in terms 
of a magnetic field configuration where rotation of magnetic field is proportional 
to the field itself thus giving rise to a constant of proportionality. The constant 
of proportionality can be a real constant (linear force free) or dependent on space 
(nonlinc;ar force free). The validity of a force free model is however restricted to 
the formation of the non evolutionary magnetic loops essentially in the vicinity 
of the solar corona. Further they are strictly invalid to describe the evolution dy- 
namics. This further necessitates the development of a self-consistent description of 
coronal magnetic field loops that can respond to as well as interact with the real- 
istic perturbations ubiquitously present in the solar atmosphere. Motivated by the 
issues described above, we have developed a more generic model based on three di- 
mensional simulations of fully compressible, non adiabatic, driven-dissipative mag- 
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netofluid plasma that contains a fairly broader spectrum of characteristic modes. 
The underlying modes range from largest scale (of the size of the system) to the 
intermediate scales (constituting the inertial range spectra) to the smallest scales 
where energy of the system can be dissipated by virtue of collisional or viscous dis- 
sipation. One of the novel features of our 3D compressible MHD plasma model is 
that it deals with the entire spectrum of fluctuations unlike those work that either 
describe single mode of fltix of coronal magnetic field and ignore the background 
small scale realistic perturbations or single coherent structures (Amari and Luciani, 
2000). Additional features of our simulation model arc that it (i) includes nonlinear 
interactions amongst a wide range of fluctuations that are initialized with random 
spectral amplitudes. The nonlinearities in the underlying system drive turbulent 
processes and lead to the cascade of spectral energy in the inertial range spectrum, 
(ii) large scale as well as small scale perturbation that may have been induced by 
the background plasma fluctuations, (iii) non adiabatic exchange of energy leading 
to the migration of energy from the energy containing modes or randomly injected 
energy driven by perturbations and further dissipated by the smaller scales. 

The plan of the paper is as follows: In section 2 we briefly summarize some of 
the MHD relaxation models appropriate for driven dissipative plasma based on 
the principle of minimum dissipation rate (MDR) leading to non-force free states. 
Section 3 presents a description of our model starting with the conservative forms of 
the MHD equations with applicable normalization procedures. Section 4 describes 
our simulation results. Section 5 include summary of our work, conclusion and 
discussions of future works, 

2. Non-force free Self-organized States from MHD Relaxation 
Models 

The seminal work on plasma relaxation proposed by Taylor predicts a force free 
state for a magnetized plasma. For such states, the current J is always along the 
direction of the magnetic field B and the perpendicular component of the current 
Jj_ = 0. So, from the force-balance equation J x B = Vp, it is seen that for force- 
free states, the pressure gradient is zero, and such states are not suitable for devices 
confining plasma by magnetic fields. 

A small amount of resistivity, ingrained in any realistic plasma, is essential to 
allow rcconncctivc processes leading to relaxation. In fact, dissipation, along with 
nonlinearity, is universal in systems evolving towards self-organized states and it 
is natural to assume that dissipation plays a decisive role in the self-organization 
of a system. Alternative models plasma relaxation, based on the principle of min- 
imum dissipation rate (MDR) of energy have been proposed by several authors 
(Montgomery and Phillips, 1988; Farengo, and Sobehart, 1994, 1995; Farengo and 
Kaputi, 2002: Dasgupta ct al, 1998, 2002, 2009; Bhattacharya and .Janaki, 2004). 
The principle of minimum dissipation rate is closely akin to the principle of min- 
imum entropy production rate of irreversible thermodynamics as formulated by 
Prigogine (Prigoginc, 1946) . Relaxation of a driven dissipative plasma has been 
formulated by Bhattacharya and Janaki (2004). As shown by these authors, the 
relaxed states obtained from MDR are non- force free, i.e., for these states, Vp ^ 0. 
A critical signature of such non-force free state is that the perpendicular component 
of the current J_l 9^ 0. It may be mentioned that relaxation of a MHD plasma to a 
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non-force free state has been numerically demonstrated by Zhu et al., (1995). One 
of the main objectives of this work is to investigate the existence of non-zero Jj_ in 
a self-organized states of a driven system. 



3. Description of the Simulation Model 

The fluid model describing nonlinear turbulent processes in the magnctofluid plasma, 
in the presence of a background magnetic field, can be cast into plasma density (pp), 
velocity (Up), magnetic field (B), pressure (Pp) components according to the con- 
servative form 

^*^'+V-Qp = Q, 



dt 



(3.1) 



where. 



F = 

■Tp 



Pp 

PpUp 
B 




6p 





I -BE 
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BU„ 



Q = 







B(Up-B) 



fM(r,i) +AtV2Up +r7V(V 
r/V^B 




Up) 



and 



7 - 1 Stt ■ 



Equations (??) are normalized by typical length £o and time to ~ f-o /Va scales in our 
simulations such that V = 4V, d/dt = tod/dt, Up = Up/Vk, B = B/yA(47^po)^/^ -P ■■ 
P/poV^, Cp = Cpl poVa,p = pI po- The bars are removed from the normalized equa- 
tions (1). Va = -Bo/(47rpo)^^^ is the Alfven speed, and I is a unit tensor. 

The rhs in the momentum equation denotes a forcing functions (fM(r, i)) that 
essentially influences the plasma momentum at the larger length scale in our sim- 
ulation model. With the help of this function, we drive energy in the large scale 
eddies to sustain the magnetized turbulent interactions. In the absence of forc- 
ing, the turbulence continues to decay freely. While the driving term modifies the 
momentum of plasma, wc conserve density (since we neglect photoionization and 
recombination). The large-scale random driving of turbulence can correspond to 
external forces or instabilities for example fast and slow streams, merged interac- 
tion region etc in the solar wind, supernova explosions, stellar winds in the ISM, 
etc. The magnetic field evolution is governed by the usual induction equation and 
obeys the frozen-in-field theorem unless nonlinear dissipative mechanism introduces 
small-scale damping. Note carefully that MHD plasma momentrum equation con- 
tains nonlinear terms on the rhs. This means that mode coupling processes can 
potentially be mediated by nonlinear interactions, in addition to the damping as- 
sociated with the small-scale turbulent motion. Thus nonlinear turbulent cascades 
are not only responsible for the spectral transfer of energy in the inertial range, but 
are also likely to damp the plasma motion in a complex manner. Nonetheless, the 
spatio-temporal scales in the nonlinear damping can be distinct from that of the 
linear dissipation. We restrict forcing of plasma momentum fluctuations in a region 
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specified by the wavenumbers k < kf such that energy is injected in the large scale 
plasma fiuctuations. The driving is random in time and space. The amplitude of 
driving force is also chosen randomly between and 1. We further make sure that 
the random number generator used for the driving force is a isotropic and uniform 
in the spectral space and does not lead to spectral anisotropy. 

Turbulent-relaxation evolution studies in three (3D) dimensions are performed 
to investigate the nonlinear mode coupling interaction of a decaying compressible 
MHD turbulence described by the closed set of Eqs (??). For this purpose, we 
have developed a full three dimensional (3D) compressible magnetohydrodynamic 
(MHD) code (Shaikh et al., 2008). All the fluctuations are initialized isotropically 
(no mean fields are asstuned) with random phases and amplitudes in Fourier space 
and evolved further by integration of Eqs (??) using a fully de-aliased pseudospec- 
tral numerical scheme. Fourier spectral methods are remarkably successful in de- 
scribing turbulent flows in a variety of plasma and hydrodynamic (i.e non mag- 
netized) fluids. Not only do they provide an accurate representation of the fluid 
fluctuations in the Fourier space, but they are also non-dissipative. Because of the 
latter, nonlinear mode coupling interactions preserve ideal rugged invariants of fluid 
flows, unlike flnite diflierence or finite volume methods. The conservation of the ideal 
invariants (energy, enstrophy, magnetic potential, helicity etc.) in turbulence is an 
extremely important feature in general, and particularly in our simulations, because 
these quantities describe the cascade of energy in the inertial regime, where turbu- 
lence is, in principle, free from large-scale forcing as well as small scale dissipation. 
The precise measurement of the decay rates associated with the MHD invariants is 
therefore one the major concerns in the study of MDR. Dissipation is nonetheless 
added physically in our simulations to push the spectral cascades further down to 
the smallest scales and also to allow minimal dissipation. The evolution variables 
are discretized in Fourier space and we use periodic boundary conditions. The ini- 
tial isotropic turbulent spectrum was chosen to be close to fc"^ with random phases 
in all three directions. The choice of such (or even a flatter than -2) spectrum do 
not influence the dynamical evolution of the turbulent fluctuations as flnal state 
in all our simulations leads to the identical results that are consistent with the 
proposed analytic theory. The equations are advanced in time using a second-order 
predictor-corrector scheme. The code is made stable by a proper de-aliasing of spu- 
rious Fourier modes and choosing a relatively small time step in the simulations. 
Additionally, the code preserves the V • B = condition at each time step. Our 
code is massively parallelized using Message Passing Interface (MPI) libraries to 
facilitate higher resolution in a 3D volume. Kinetic and magnetic energies are also 
equi-partitioned between the initial velocity and the magnetic fields. The latter 
helps treat the transverse or shear Alfven and the fast / slow magnetosonic waves on 
an equal footing, at least during the early phase of the simulations. 



4. Simulation Results: 

Our major focus in the paper is to study turbulent relaxation of magnetofluid 
plasma through nonlinear interactions. For this purpose, we let magnetized fully 
compressible MHD turbulence evolve under the action of nonlinear interactions in 
which random initial fluctuations, containing sources of free energy, lead to the 
excitation of unstable modes. These modes often deviate the initial evolutionary 
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system substantially away from its equilibrium state. When the instability tends to 
saturate, unstable modes lead to fully developed turbulence in which larger eddies 
transfer their energy to smaller ones through a forward cascade until the process is 
terminated by the small-scale dissipation. During this process, MHD turbulent fluc- 
tuations are dissipated gradually due to the finite Reynolds number, thereby damp- 
ing small scale motion as well. The energy in the smaller Fourier modes migrates 
towards the higher Fourier modes following essentially the vector triad interactions 
k + p = q. These interactions involve the neighboring Fourier components (k, p, q) 
that are excited in the local inertial range turbulence. We conjecture in our earlier 
work (Shaikh et al., 2008) that MDR plasma relaxes towards nonlinear non force 
free state state through nonlinear evolution. In our simulations, we undertake this 
point in the presence of driven turbulence. Our simulations arc carried out in the 
presence of a mean magnetic field (taken along the z direction). One of the ways we 
ensure the nonlinear non force free evolution is by determining whether or not there 
develops any perpendicular component of plasma currents. In the following, we ex- 
plain why it is essential to self-consistently produce a perpendicular component of 
plasma current ( J_l) that facilitates the nonlinear non force free interactions. 

The nonlinear interactions in magnetoplasma turbulence arc; detcirmincid pre- 
dominantly hy J X B forces in plasma momentum equation. Similarly, since V oc B 
through (J), the same nonlinear interactions also govern magnetic induction equa- 
tion. In the presence of a mean magnetic field along the z direction, the parallel 
and perpendicular components can be denoted respectively by J|| and J±. For ob- 
vious reasons, J|| x B = for magnetic field fluctuations that lie strictly along the 
the J|| component of plasma current. Hence this term contributes negligibly in the 
nonlinear interactions. By contrast, it is only the J± component, emerging from 
the J± X B term, that contributes largely to the nonlinear interactions. Thus cur- 
rent fluctuations orthogonal to the mean or fluctuating magnetic field component 
predominantly govern the nonlinear interactions. 

Figure ?? shows iso-surfaces of the a;-component of the magnetic field by the 
evolution of random initial turbulent fluctuations leading to the formation of rel- 
atively small-scale isotropic structures. Fig 1 is a typical snap shot of nonlinear 
turbulent fluctuations during an early phase of evolution. The initial condition in 
combination with the random forcing leads to the state depicted as figure 1. In 
Figure ??, we follow the evolution of both Jy and J± components of the currents to 
quantitatively measure their progressive development. Clearly, the two components 
gradually develop and evolve self-consistently in our simulations. This further leads 
to Jj_ = J — (J • B/|B|)6 - where b = B/|B| is the unit vector along B, - which 
is non-zero, thus proving conclusively that the resultant state is non force-free. For 
force-free state, J± is zero 

The corresponding decay rates associated with turbulent relaxation of the rugged 
ideal invariant of MHD, viz, magnetic helicity {K = J A • B dv), magnetic energy 
{E = 1/2 J fiv) and energy dissipation ratc(7? = rj J rfv) are shown simulta- 
neously in Figure ??. Clearly, the magnetic energy E decays faster than magnetic 
helicity K, and energy dissipation (R) decays even faster than the two invariants. 
This state corresponds to a minimun dissipation in which selective decay processes 
lead to the faster decay rates of the magnetic energy (when compared with the 
magnetic helicity decay rates). Furthermore, the time evolutions of the volume av- 
erages of global heliciy, magnetic energy and the dissipation rates are plotted in 
Figs 3. Thus, these quantities are averaged over the fluctuations, and hence they 
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show a regular behavior. So the decay of the global heliciy, magnetic energy and the 
dissipation rates are not due to any linear decay they are the results of nonlinear 
turbulence in the system. The selective decay processes (in addition to dissipations) 
depend critically on the cascade properties associated with the rugged MHD invari- 
ants that eventually govern the spectral transfer in the inertial range. This can be 
elucidated as follows (Biscamp, 2003). Magnetic vector potential in 3D MHD dom- 
inates, over the magnetic field fluctuations, at the smaller Fourier modes, which 
in turn leads to a domination of the magnetic helicity invariant over the magnetic 
energy. On the other hand, dissipation occurs predominantly at the higher Fourier 
modes which give rise to a rapid damping of the energy dissipative quantity R. A 
heuristic argument for this process can be formulated in the following way. The 
decay rates of helicity K and dissipation rate R = Jy r/j^dV in the dimcnsionless 
form, with the magnetic field Fourier decomposed as B(k, i) = I]fcbkexp(ik • r), 

dK 27J 2 dR 2if . 2 

— =--E.fcbk; _ = -_E.fcbk (4.1) 

where S = T^|TA^ is the Landquist number, tjj,ta arc the resistive and Alfven 
time scales, respectively. The Landqiiist mimbor in our simulations varies between 
10^ and 10^. We find that at scale lengths for which k w S*^/^, the decay rate 
of energy dissipation is ^ 0(1) 1 . But at these scale lengths, helicity dissipation 
is only ~ 0{S~^^^) ^ 1. This physical scenario is further consistent with our 3D 
simulations. Interestingly, the decay rates of kinetic energy of turbulent fluctuations 
are initially higher than the magnetic energy. Hence the ratio of magnetic to kinetic 
energies shows a sharp rise in the initial evolution, as shown in Figure ??. However, 
as the evolution progresses, the two decay rates become identical and the ratio 
approaches eventually a constant value. Another important outcome to emerge from 
our investigations is that a state corresponding to the minimum dissipation rates, is 
more plausible in a driven dissipative plasma. So, we may conclude notwithstanding 
the Taylor hypothesis of force-free state, our simulations clearly demonstrate that 
the nonlinear selective decay processes lead the plasma fluctuations to relax towards 
a non-force-free state in a rather natural and self-consistent manner. 

5. Summary, Conclusion and Future Work 

We have demonstrated through a fully self-consistent 3-D numerical simulation that 

a compressible, dissipative magneto-plasma driven by large-scale perturbations, can 
give rise to a self-organized state, which is not force-free. This is corroborated by 
the presence of a perpendicular component of the current J± in the self-organized 
state. In addition, the comparative decay rates of global helicity, magnetic energy 
and the (Ohmic) dissipation rate amply indicate that dissipation rate can also serve 
as an effective minimizer for a driven dissipative plasma. 

It is to be noted that the imderlying system has a finite dissipation. When the 
rate of dissipation exceeds that of the forcing, dissipative processes dominate the 
evolution. In such a case, both the energy and helicity decreases rapidly. Note that 
dissipation is necessary in our simulation to validate the hypothesis of minimum 
dissipation rates that lead eventually to a self-organization state in a dissipative 
MHD plasma by selectively operating on magnetic energy and helicity. In the event 
of initially zero fluctuations, the forcing leads to populate the turbulent spectra 
during the early phase of evolution, which then decay owing to the flnite dissipation 
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in the system. The final results in such cases give rise to the same conclusions that 
are described in the context of MDR state elsewhere in our paper. 

Further, we would like to indicate that our results are in line with Prigoginc's idea 
of the emergence of new type of self-organized states called 'Dissipative Structures' 
for systems far from equilibrium, where the non-linear interactions of fluctuations 
create ordered state in the system. A detailed investigations on the dissipative 
structures in plasma will be undertaken as our future works. 
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